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ABSTRACT 

It is believed that the local massive black holes were dominantly grown up through accretion during 
quasar phases, while a fraction of the local black hole mass was accumulated through accreting gases 
at very low rates. We derive the black hole mass density as a function of redshift with the bolometric 
luminosity function of active galactic nuclei (AGN) assuming that massive black holes grew via ac- 
creting the circumnuclear gases, in which the derived black hole mass density is required to match the 
measured local black hole mass density at z = 0. Advection dominated accretion flows (ADAFs) are 
supposed to present in low luminosity active galactic nuclei (AGNs)/normal galaxies, which are very 
hot and radiate mostly in the hard X-ray band. Most of the X-ray background (XRB) is contributed 
by bright AGNs, and a variety of AGN population synthesis models were developed to model the 
observed XRB in the last two decades. Based on our derived black hole mass density, we calculate 
the contribution to the XRB from the ADAFs in faint AGNs/normal galaxies with a given Eddington 
ratio distribution, which is mostly in hard X-ray energy band with an energy peak at r^ 200 keV. The 
growth of massive black holes during ADAF phase can therefore be constrained with the observed 
XRB. Combining an AGN population synthesis model with our results, we find that the fitting on 
the observed XRB, especially at hard X-ray energy band with > 100 keV, is improved provided the 
contribution of the ADAFs in low-luminosity AGNs/normal galaxies is properly included. It is found 
that less than ~15 per cent of local massive black hole mass density was accreted during ADAF phases. 
We suggest that more accurate measurements of the XRB in the energy band with > 100 keV in the 
future may help constrain the growth of massive black holes at their late stage. We also calculate 
their contribution to the extragalactic 7-ray background (EGRB), and find that less than ^ 1% of 
the observed EGRB is contributed by the ADAFs in these faint sources. 

Subject headings: galaxies: active — quasars: general — accretion, accretion disks — black hole physics; 
X-rays: diffuse background 



1. INTRODUCTION 

It is well believed that almost all galaxies contain 
massive black holes at their centers, and a tight cor- 
relation was revealed between central massive black 
hole mass and the veloc i ty dispersion of the galaxy 
(jFerrarese fc Merritli [2OOOI : iGebhardt et al.|[2000l ). The 
growth of massive black holes at the centers of g alaxies 
may p robably be linke d to a ccretion processes (jSoltanl 
II 9821 ). lYu fc Tremaind ()2002D estimated the black hole 
masses from the stellar velocity dispersions of galaxies 
measured by the Sloan Digital Sky Survey (SDSS) us- 
ing the empirical relation between black hole mass and 
the velocity dispersion, and the local black hole mass 
density was derived. They further calculated the black 
hole mass density accreted during optical bright quasar 
phases using an optical quasar luminosity function (LF), 
and found that the accreted mass density is consistent 
with the local black hole mass density estimated from 
the velocity dispersions, if a radi ative efficiency ^ 0. 1 
is adopted for quasars (also see [M arconi et al.' 12004 



Shankar et al]l2004t iHopkins et alt l2007: Shanka r et al 



2009( 1. This implies that the growth of massive black 



holes through accretion during optically bright quasar 
phases may probably be important, if a radiative effi- 
ciency ^ 0.1 is adopted. If the massive black holes are 
spinning rapidly, their radiative efficiency can be higher 



than that they adopted, the black hole mass density ac- 
creted during quasar phases would be lowe r than the 
measured local black hole mas s density fe.g.. lCaoll2007l : 
iCao fc Lil [200a iLiet al.l[2010l ). Thus, one cannot ne- 
glect the contribution from accretion at low rates to the 
growth of massive black holes, provided the duration of 
accretion at low rates is as long as Hubble timescale. 

The cosmological X- ray background (XRB) is mostly 
contr ibuted by AGNs (|Hasingeii Il998t iSreekumar et al.1 
119981 ). which can be used t o constrain mass ive black 
hole accretion history (e.g.. lElvis et al.ll2002| ). In the 
most popular synthesis models of the XRB based on the 
unification schemes for AGNs, the cosmological XRB 
contributed by Compton-thin AGNs can account for 
~ 80% of the observed XRB (e.g., lUeda et all [20031: 
[Gilhet al.l [20071) . The residual XRB can be explained 
provided the same number of Compton-thick AGNs with 
log N-a = 24-25 as thos e with log JVh = 23 - 2 4 is iii - 
cluded ([Uedaet al.ll2003f ). Idi Matteo fc FabianI (|1997[ ) 
have alternatively proposed that the hard XRB above 10 
keV may be dominated by the thermal bremsstrahlung 
emission from the advection dominated accretion flows 
(ADAFs) in low-luminosity AGNs. Due to the diffi- 
culties on detecting Compton-thick AGNs, the space 
number densit y of Compton-th ick AGNs is still quite 
debated (e.g., lUeda et al.l [2003t [Treister fc Ifa^ [200I 
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iGilli et al.|[2007l: iTreister et al.|[2009l ). Bas ed on the sur- 
veys w ith Swift and INTEGRAL satelhtes. lTreister et alj 
(|2009f ) suggested that the space number density of 
Compton-thick AGNs may be significantly lower than 
those adopted in the previous AGN population synthe- 
sis models to explain the observed XRB. Due to the 
uncertainty of the number density of Compton-thick 
AGNs, the residual XRB could be attributed to both 
the Compton-thick AGNs and the ADAFs in faint AGNs 
and normal galaxies. Therefore, the observed XRB can 
be used to constrain massive black hole accretion history 
at their late stag e when they are accreting at low rates 
(|Caol 120051 [2007fl . ICaol (l2007l ) calculated the contribu- 
tion to the XRB from the ADAFs in faint AGNs/normal 
galaxies, and compared it with the residual XRB. Their 
results showed that less than ~5 per cent of local mas- 
sive black hole mass density was accreted during ADAF 
phases, otherwise the XRB contributed from ADAFs will 
surpass the observed residual XRB even if no Compton- 
thick AGNs are included. For simplicity, they adopted 
an average mass accretion rate for faint AGNs/normal 
galaxies in calculating the contribution of ADAFs to the 
XRB. The radiative efficiency of standard thin accre- 
tion disks in bright AGNs does not vary with mass ac- 
cretion rate. However, the faint AGNs/normal galaxies 
may probably contain ADAFs, of which t he radiative ef- 
ficien cies vary with mass accretion rate (jNaravan fc Yil 
11995( 1. The average radiative efficiency for a popula- 
tion of sources containing ADAFs accreting at different 
rates can be calculated by weighing over the distribution 
of mass accretion rate. The Eddington ratios of AGNs 
spread over several orders of niagnitu de fe.g.. iHol 120021 : 
IHopkins et all [2006t ICao fc Xul 120071) . The Eddington 
ratio distribution for accreting massive black hole holes 
derived from observations exhibits nearly a power-law 
distri bution with an exponential cutoff at high Eddington 



xpc 
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ratio (iMerloni fc Heinz)l2008l: IHopkins fc Herna'^Iistl l2009l: 
iKauffmann fc: HeckmanI l2009f K In this paper, we will 
ad opt the Edd ington ratio distribution of AGNs derived 
bv IHopkins fc Hcrnauist (2009) to calculate the contri- 
bution of faint AGNs/normal galaxies to the XRB. Com- 
pared with the residual XRB, the constraints on the frac- 
tion of local black hole mass accreted in ADAF phases 
are derived. The cosmological parameters flu = 0.3, 
n\ — 0.7, and Hq = 70 km s~^ Mpc~^ have been 
adopted in this work. 

2. BLACK HOLE MASS FUNCTIONS 

In this section, we derive the black hole mass den- 
sities as func t ions o f redshift in a similar way as that 
done in iCaol (|2007| ). In this work, we use the bolo- 
metric quasar lumin osity function (QLF) derived by 
IHopkins et al.l (|2007f ) to calculate the black hole mass 
densities. We summarize our calculations as follows (see 
[Cao 2007, for the details). 

iHopkins et al. (2007 ) 's QLF is calculated by using a 
large set of observed quasar luminosity functions in var- 
ious wavebands, from the IR throug h optical, soft and 
hard X-rays (see IHopkins et al1l2007l for the details). 
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with normalization (/>*, break luminosity _L», faint-end 
slope 7i, and bright-end slope 72. The break luminosity 



L, evolves with the redshift is given by 

logi, = (logL,)o + kL,l^ + /CL,2^^ + fcL,35^ (2) 

and the two slopes 71 and 72 evolves with redshift as 

71 = (7i)o(- ^^''^1- 
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The parameter ^ is 



e = iog( 
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and Zj-ef = 2 is fixed. The best-fit parameters we adopted 
in this work are as follows: log(?!)*(Mpc~^) = —4.825 ± 
0.060, [logL43.9 X 1033ergs-i)]o = 13.036 ± 0.043, 
fcii = 0.632 ± 0.077, kL2 = -11.76 ± 0.38, /cl 3 = 
-14.25 ±0.80, (71)0 = 0.417 ±0.055, kj, = -0.623 ± 
0.132, (72)0 = 2.174 ± 0.055 , fc^,.i = 1.460 ± 0.096, and 
fc^^_2 = -0.793 ± 0.057 ( Ho pkins et al.|[2007l ). 

In this work, all the sources described by this QLF 
$(L, z) are referred as active galaxies, which are sup- 
posed to have mass accretion rates to > TOcrit = 0.01 
[defined as m = O.lMc^/iEdd, and LEdd — 1-3 x 
2o38 g].g s~-'^(Mbh/M0)] and contain standard radiative 
efficient accretion disks. The cosmological evolution of 
black hole mass density caused by accretion during ac- 
tive galaxy phases is described by 



dpt^jz) _ dt 1 
dz dz M0 



(1 -e)Ld$(L,z) 
ec^ d log L 



dlogi, (6) 



where ^{L, z) is the bolometric QLF given by 
IHopkins et all (|2007t ). p^^^{z) (in units of Mq Mpc'^) is 
the black hole mass density accreted during active galaxy 
phases from Zmax to z, and e is the average radiative ef- 
ficiency for active galaxies. 

For those inactive galaxies with L < 10"*^ ergs s~^, 
their mass accretion rates are very low. The term "inac- 
tive galaxies" used in this work does not mean that they 
are really inactive. Compared with active galaxies, the 
massive black holes in inactive galaxies are still accreting, 
but at low rates with m < TOcrit = 0.01, and radiatively 
inefficient ADAFs are suggested to be present in inac- 
tive galaxies (e.g.. lNaravanll2002l ). The black hole mass 
density Pbh(z) accreted during inactive galaxy phases be- 
tween z and Zmax can be calculated by 



dp^^_pir*(fK 
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where TOf™^^ is the average dimensionless mass accretion 
rate for the inactive galaxies, pjX'^'(z) is the black hole 
mass density of inactive galaxies at redshift z, g^DAF jg 
the average radiative efficiency of the ADAFs in those 
inactive galaxies. The radiative efficiency e^^AF jg ^g^_ 
ally much lower than 0.01, so we approximately adopt 
gADAF ^ Q JQ Q^j. calculations on black hole mass densi- 
ties. 

Assuming that the growth of massive black hole is dom- 
inated by accretion, the total black hole mass density is 
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given by 



Pbh(^)+Pbh(^)+Pbh(^;max), 

(8) 



where Phhizmux) is the total black hole mass density at 
Zmax- For active galaxies, the black hole mass density at 
redshift z can be calculated from the QLF by 



exponential term in (jl2[) . which is a good approximation 
because m < rhcrit << Vi ^^'^ the accretion rate distribu- 
tion (J12p is always normalized by assuming all inactive 
black holes to be accreting with rates in the range of 



< m < TOn 



Thus, the accretion rate 



Pbh(^) 
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act ^ 



Edd,0 



^^^^dlogL MoMpc-3, 
dlogL 



(9) 

where $(L,z) is the bolometric QLF given by Eq. ([T]), 
and A^ct' — -^boi/^Edd is the average Eddington ra- 
tio for active galaxies. In this work, we adopt A^c^ = 
0.25 as that derived from a large bright AGN sample 
(jKoUmeier et al.l 120061 ) . The black hole mass density for 
inactive galaxies p^^'^^iz) can be calculated with 



distribution can be described by one parameter Km with 
specified mass accretion rate range for inactive galax- 
ies. The standard thin accretion disks are present in 
bright quasars, while they will transit to ADAFs pro- 
vided the mass accretion rates are lower than the critical 
value mciit (e.g ., Narayan. Mahadcvan & Quataert 199g; 
lNaravanll2002f ). The critical dimensionless mass accre- 
tion rate merit — 0.01 is suggested eit her by observat ions 
or theoretical model calculations (see lNaravanll200"^ for 
a review and references therein). For inactive galaxies 
with m < TOcrit, their average mass accretion rate can be 
calculated with 



pkr'w = pbhW-<(^) 



Pbh(2)+Pbh(2)+Pbh(Zmax)- 
(10) 

We assume pl^iz) = 0.5/9bh(^max) at z = z^ax, i-e., 
a half of all massive black h oles a re active at Zrnax, as 
that used in iMarconi et al.l (J200J) , and the black hole 
mass density Phh{zmax) = '^Phhi^ma.x) is obtained from 
Eq. ([9]). Thus, the black hole mass density as func- 
tion of redshift z can finally be calculated by integrat- 
ing Eqs. ([6]) and (O from Zmax to z with Eqs. ([8])- 

and e are 



act 
Pbh 



(z). 
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act 
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specified. We adopt Zmax = 5 in all calculations, and 
the value of e is tuned to make the derived total black 
hole mass density phh{z) match the local black hole mass 
density p]^^''Hz) = 4.0 x lO^M,^ ^V^' ^ at z = (e.g., 
iSalucci et aIlll999HMarconi et al.ll2004[) . when the value 
of TOf™^^. is specified. The black hole mass accreted dur- 
ing z > 5, Phhiz = 5), should be ^ /5bii(0), which implies 
that it will make little difference on our calculations even 
if a better QLF is available for Zmax > 5. Thus, we will 
not extrapolate the present QLF to higher redshifts. 

3. EDDINGTON RATIO DISTRIBUTION FOR INACTIVE 
GALAXIES 

The observed Eddington ratio distribution can be de- 
scribed by 



/a (A) 



diV 



A^dlogA 



Co 



exp 



(11) 



where A is the Eddi ngton ratio, A = L/ Ly.ru, and Cq 
is the normalization (jHopkins fc Hernauistll2009() . They 
suggested that k w 0.3 - 0.8, and ?? = 0.2 - 0.4. This 
Eddington ratio distribution is consistent with the self- 
regulated black hole growth model, in which feedback 
produces a self-regulating "decay" or "blowout" phase 
after the AGN reaches some peak lumi nosity and be- 
gins to e xpel gas and shut down accr etion (jHopkins et ahl 
I2005allbl : iHopkins fc Hernaulstl[200l . 

In this work, we assume the distribution of dimension- 
less mass accretion rate have a similar form as Eq. ([TT|) . 



/m(TO) = 



diV 



A^d log m 



CA- 

V 



exp 






(12) 

for low-luminosity AGNs with m < merit, where C2 is 
the normalization. In all our calculations, we drop the 



''inact 



m/m(m)dlogm, 



(13) 



where the minimum accretion rate, rrtnun = 1.0 x 10~^ 
is adopted in all our calculations. Thus, the value of 
parameter k^ corresponds to an average mass accretion 



rate m-- . 



4. CONTRIBUTIONS OF INACTIVE GALAXIES TO THE 
COSMOLOGICAL BACKGROUND RADIATION 



We employ the approach suggested by iManmotol 
()2000[) to calculate the global structure of an ADAF sur- 
rounding a massive black hole in the general relativis- 
tic frame. All the radiation proces ses are included in 
the global structure calculations (see lManmotoll2000l for 
details and the references therein). The global struc- 
ture of an ADAF surrounding a IO^Mq black hole with 
spin parameter a can be calculated, if the model pa- 
rameters, dimensionless mass accretion rate m, mag- 
netic field strength relative to gas pressure /3, defined as 
Pm = -B^/Btt = (1 - l3)ptot [ptot = Pgas +Pm), the frac- 
tion of the released gravitational energy directly heating 
the electrons 6, and the conventional viscosity parameter 
a, are specified. The constr aints on the values of these 
parameters were discussed in lCaol (|2007D . We adopt the 
same values as that work, i.e., a = 0.2, and /3 = 0.8, 
with which no global solution is available for m > 0.01. 
This is consistent with the observations (e.g.. iNaravanI 
l2002f ). The value of S is still a controversial issue. In 
most of our calculations, we adopt a conventional value 
of 5 = 0.1 as that adopted ICaol (|2007D . We also calculate 
the case with 6 = 0.01 for comparison. 

The inner region of the ADAF is very hot, and the tem- 
perature of protons can be as high as ^ lO^^K. Thus, 
7-ray emission may be produced through the pion pro- 
duction processes in the proton-proton (p-p) collisions 
and subseque ntly decay of neutral p ions, which can be 
described by (jMahadevan et al.|[T997[ ) 



P + P 



P+P + TT°, TT° 



71 +72- 



(14) 



The production of gamma rays thro ugh p-p colli- 
sions has been studied in s everal works ( Dermeii [T986I : 
iGiovanneUi et al.lll982al lH). iSteckeii U22M showed that 
the gamma-ray spectrum produced in unit volume of the 
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flow from 7r° decay can be calculated by 



L^ = h^v f^{E^) ^2h^v 



dK 






m^ 



(15) 

where E^^ is the pion energy in GeV, TTi^r is the mass of 
the pion in GeV c~^, /ttI-Ett) is the tt" spectrum, and 
^TTmin ^^ thc minimum pion energy required to produce a 
gamma-ray with energy E^ , which is described by 



-E-TT 



-E-, 



m: 
4E. 



(16) 
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The gamma-ray spectrum of an ADAF mainly depends 
on the density and temperatu re distribution s of the ions. 
The n^ spectrum is given by ()Dermei1ll986[ ). 



/ttI-Ett) 



cn„ 



d7r 



(7.^ - 1) 



im^9pKl{l/9p) J, "^ [2(7, + 1)] 1/2 



X ;'^""d7*(r7*)-^^"*^^*'^'-^ 



d7* 



x{exp[-977*(l-/3/3*)]-exp[-q77*(l+/3r)]} 
cm-3 s-i GeV~\ (17) 

where Up is the number density of ions, 0p = kTi/iripC^ is 
the dimensionless ion temperature, E^^ = ^niTr is the pion 
energy in the observer's frame, and K2{x) is the modified 
Bessel function of order 2. The differential cross section 
for the production of a neutral pion with Lorentz factor 
7* in the center-of-momentum system (CM) of two col- 
liding protons with relative Lorentz factor 7r is de noted 
by dcr*(7*; 7r)/d7*, which is given in lSteckeii()1971[ ). The 
quantities q and C, are defined as q = [2(7 -I- l)]i/^/6'p 
and C = {S - Ami + ml)/2S^I'^, where S = 2m2(7, + 1) 
(see lSteckei1ll97ir for the details). Thus, both the X-ray 
and 7-ray spectra of an ADAF can be calculated when 
the global structure of the ADAF is available. 

We perform a set of spectral calculations for ADAFs 
surrounding black holes with different masses, and find 
that the dependence of ADAF spectrum on black hole 
mass is almost perfectly linear in X-ray/7-ray bands for 
-^^bh ^ IO^Mq, provided all other parameters are fixed. 
Therefore we simply use the spectrum lEirn) of an ADAF 
around a typical massive black hole, lO^M©, accreting at 
the rate m as a template spectrum to calculate the con- 
tribution of inactive galaxies in unit of co-moving volume 
to the XRB by muhiplying pj^'jf'=*(z)/10^MQ with the ac- 
cretion rate distribution (fT2|) (see Eq. [19] in this section). 
We can calculate the average spectrum of a population 
of ADAFs accreting at rates with a distribution given by 

(USD, 



Le{E) 



fm{m)lE{m)dlogm, 



(18) 



where Isifn) is X-ray and 7-ray spectrum from an ADAF 
surrounding a lO^M© black hole accreting with m. 

The contribution of the ADAFs in all inactive galaxies 
to the cosmological background radiation can be calcu- 



lated by 



GeV S-' Hz 



'^ cm ^, 



fiE) 



lO^M 



Jo 



Pwr*w(i + ^)^E[(i + ^)E]di/ 



4ndl 



dz, 



dz 
(19) 

where Le{E) is the template spectrum of ADAFs sur- 
rounding 1O^M0 black holes averaged over accretion rate 
rh. 



5. RESULTS 

We plot average dimensionless mass accretion rate 
"^fnact ^ ^ function of power Km of the accretion rate 
distribution (see Eq. [T2l) for the inactive galaxies in Fig. 
[TJ As discussed in §2, we can calculate both active and 
inactive black hole mass densities as functions of red- 
shift simultaneously using the bolometric QLF based on 
the assumption that the growth of massive black holes 
is dominated by mass accretion. We plot the total black 
hole mass density pbh(2), the black hole mass densities 
Phh^^i^) foi' inactive galaxies and p^h_{^) ^^^ active galax- 
ies, as functions of redshift z in Fig. [21 for different values 
of Km (i.e., to|™''), respectively. The ratios of black hole 
mass densities P^y,{z) accumulated during ADAF phases 
between z and Zmax to the local black hole mass density 
^locai ^ 4 X IO^MqMpc"^ are also plotted in Fig. [1 
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Fig. 1. — The relation of the power Km in Eq. I I12I I with average 
dimensionless mass accretion rate ni?™^^ for inactive galaxies. The 
minimum accretion rate, rhmin = 1-0 x 10~^ i s ad opted in the 
calculations. The relation of the power Km in Eq. I I12II with average 
dimensionless mass accretion rate 'rnf^^'^^ for inactive galaxies. The 
minimum accretion rate, rhmin = 1-0 X 10~^ is adopted in the 
calculations. 
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Fig. 2. — The total black hole mass densities Pbh(^) accumulated 
through accretion as functions of redshift z (solid lines) in unit of 
^local _ 4 X IO^MqMpc"^. The dashed and dotted lines represent 
the black hole mass densities p^^^'^^(,z) for inactive galaxies and 
pj^^'(z) for active galaxies at redshift z, respectively. The dot- 
dashed lines represent the ratio of black hole mass densities (^^ (z) 
accumulated during ADAF phases to the local black hole mass 
density Pbij"""' = 4 X IO^MqMpc"^ from Zmax to z. The different 
colors correspond to different values of Km = 0.2 (black) and 0.3 
(red). 

The global structure of an ADAF surrounding a spin- 
ning massive black hole is available by solving a set 
of genera l relat ivistic hydrodynamical equations (see 
iManmotol I200CII . for the details). The structures of 
ADAFs are plotted in Fig. [3] for different parameters, 
which show that the proton temperature in the inner 
edge of the ADAFs surrounding rapidly spinning black 
holes can be as high as ^ lO^^K. The spectra of ADAFs 
surrounding massive black holes can be calculated with 
the derived global structure of ADAFs (see §4) . In Fig. 
m we plot the X-ray and 7-ray spectra of the ADAFs 
surrounding massive black holes with different values of 
spin parameter a. We also plot the spectra calculated 
with (5 = 0.01 in Fig. [5j In the calculations, the black 
hole mass Mbh = 1O^M0 is adopted. As the radiative 
efficiency of an ADAF is no longer constant, we calcu- 
late the bolometric luminosities of ADAFs as functions 
of mass accretion rate m. Figure [5] shows how the Ed- 
dington ratios iboi/-^Edd vary with mass accretion rate 
m for different black hole spin parameters and ADAF 
parameters. The average X-ray and 7-ray spectra from a 
population of ADAFs, of which the mass accretion rate 
distribution is described by Eq. (fT2|) . are given in Fig. [7] 
for (5 = 0.1 and 0.01, respectively. 

The average mass accretion rate mf™ct can be calcu- 
lated with Eq. (J13p for a specified index Km of the power 
law accretion rate distribution (fT2|) . For a given aver- 
age mass accretion rate JTifnacf t^^ contribution of the 
ADAFs in all inactive galaxies to the X-ray /7-ray back- 
ground is calculated with the derived inactive black hole 
mass density and ADAF spectra. We plot the contribu- 
tion of the ADAFs in all inactive galaxies to the X-ray 
band cosmological background with different black hole 
and ADAF parameters in Figs. |8] and [9l For compar- 
ison, we also plot the observed XRB in the figure, and 
the sum of the contributions from the type I/II bright 



AGNs(Compton-thin) given bv lTreister et all ()2009[ ) and 
all inactive galaxies derived in this work. Their contri- 
bution to the extragalactic 7-ray background is plotted 
in Figs. [TOl and fTTJ for (5 = 0.1 and 0.01, respectively. 
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Fig. 3. — The structure profiles of the ADAFs surrounding a 
10*A/q black hole with different spins, a = 0. (blue lines), a = 
0.5(green lines) and a = 0.9(red lines) accreting at m = 0.001. 
Left panel: the temperature distributions of the electrons (dotted 
lines) and ions (solid lines) as functions of dimensionless radius 
R/Rb (Rs = 2GMI(?). Right panel: the radial velocity (dotted 
lines) and sound speed (solid lines) distributions. The parameter 
5 = 0.01 is adopted in the calculations. 




Fig. 4. — The X-ray /7-ray spectra of the ADAFs surrounding a 
IO^A/q Schwarzschild (blue lines) or Kerr black hole with a = 0.5 
(green lines) and a = 0.9 (red lines) accreting at different rates. 
The solid, dashed, and dash-dotted lines are corresponding to the 
accretion rates m = 1.0 X 10~^, 1.0 X 10""*, and 1.0 X 10^^, re- 
spectively. The parameter S = 0.1 is adopted in the calculations. 
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Fig. 5. — The same as Fig. [J] but the parameter 5 = 0.01 is 
adopted. 
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Fig. 6. — The relation between dimensionless mass accretion rate 
rh and Eddington ratio Lbol/^Edd for ADAFs in inactive galax- 
ies. The dotted and dashed lines represent the results for <5 = 0.1 
and 0.01, respectively. The lines with different colors correspond 
to different values of a = (blue), 0.5 (green), and 0.9 (red), re- 
spectively. For comparison, we also plot several power-law lines 
^bol/^Bdd oc rh" in the figure. The black dash-dotted line repre- 
sents s = 1, i.e., constant radiative efficiency, and e = 0.1. 
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Fig. 7. — The X-ray and 7-ray spectra averaged over a population 
of ADAFs surrounding IO^Mq black holes, which has rh distribu- 
tion as indicated in Eq. 1121 with different values of spin parameter 
a and ADAF model parameter S. The different color lines repre- 
sent the cases with spin parameter a = (blue), 0.5 (green), and 
0.9 (red), respectively. The solid and dashed lines correspond to 
the spectra with S = 0.1 and 5 = 0.01, respectively. The parameter 
Km = 0.3 is adopted in the calculations. 
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Fig. 8. — The contribution of the ADAFs in all inactive galax- 
ies to the XRB. The thick black solid line represents the XRB 
spectrum from the AGN population synthesis model given by 
ITreister et al.l IJ2009I ) . The thin solid line represents the contri- 
bution to their model from Compton-thick AGNs. The lower blue, 
green and red lines correspond to the contribution from all inac- 
tive galaxies with spin parameter a = 0, 0.5, and 0.9, respectively. 
The upper blue, green, and red lines correspond to the sum of X- 
ray cont ributions from both the AGN population synthesis model 
given bv lTreister et all f2009, ) and all inactive galaxies. The dashed 
and dotted lines represent different values of Km: 0.2 and 0.3, 
which correspond to average mass accretion rate of inactive galaxies 
rhf^lL- = 8.35 x 10"" and 6.12 X 10"*, respectively The param- 
eter = 0.1 is adopted in the calculat ions. The measur ements of 
the extr agalactic X RB with Chandra ( Hickox Sz Markcvitch 200(3}, 
XMM IIDe Luca fc Molcndi 2004), INTEGRAL (Churazov et at] 
1200 7) and Swift l(Siello et al.l2008l ) are plotted as blue shaded area, 
yellow shaded area, red and green points, respectively. 
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Fig. 9. — The same as Fig. \8\ but the parameter <5 = 0.01 is 
adopted. 
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Fig. 11. — The same as Fig. 1101 but the parameter <5 = 0.01 is 
adopted. 
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Fig. 10. — The contribution of the ADAFs in all inactive galax- 
ies to the EGRB. The black solid line denotes the observed EGRB 
(taken from Sreckumar et al. 1998). The different color lines corre- 
spond to the contribution from ADAFs surrounding black holes in 
all inactive galaxies with spin parameter a = (blue), 0.5 (green), 
and 0.9 (red), respectively. 



6. DISCUSSION 

Unlike standard thin accretion disks, the radiative ef- 
ficie ncy of ADAFs incre ases with mass accretion rate 
m (jNaravan fc Yil 119951) . Thus, the average radia- 
tive efficiency for a population of ADAFs is avail- 
able only if the distribution of mass accretion rate 
m is known. The Eddington ratio distribution for 
AGN s was extensively explored by different authors 
(e.g..lMerloni fc Heinzl'l200l iHopkins fc Hernauistl [2009l: 
iKauffmann fc HeckmanI 120091 )^ These works suggested 
that more sources are accreting at lower rates, and the 
Eddington ratios for AG Ns can be well described b y 
a power law distribution (jHopkins fc HernquistI 120091 ) . 
This implies that the average radiative efficiency for the 
ADAFs with this Eddington ratio distribution can be 
lower than that for an ADAF accreting at the average 
mass rate, because the radiative ef ficiency of ADAFs in - 
crea s es wit h mass accretion rate m (jNarayan fc Yil[l995| ) . 
iCaol ()2007[ ) calculated the contribution of ADAFs to the 
XRB, and the constraint on the fraction of black hole 
mass accreted in ADAF phases was derived by compar- 
ing with the resid ual XRB. A n average mass accretion 
rate is adopted in iCaol (|2007[ )'s calculations for simplic- 
ity, and therefore the upper limit on the mass growth 
of black holes in ADAF phases is under -estimated. In 
this work, we improve the calculations by ICaol (|2007D by 
including the Eddington ratio distribution for AGNs in 
our calculations. 

Using the method described in §2, the active, inac- 
tive and total black hole mass densities as functions of 
redshift are calculated by tuning the average radiative 
efficiency e for active galaxies to let the total black mass 
density match the measured local black hole mass den- 
sity at z = (see Fig. [2]). It is found that most of the 
local black hole mass was accumulated during z < 2, 
and fraction of the local black hole mass accreted during 
ADAF phases is determined by the average mass accre- 
tion rate mf^J^., i.e., the distribution of accretion rates 
for inactive galaxies. We find that the average radiative 
efficiency ci 0.134 is required for active galaxies, which 
corresponds to a « 0.9. This is consistent with that de- 
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rived in|Ci3 (|2007l ) and lElvis etall (|200^ . 

In the calculation of the black hole mass densities, 
we need to know the value of the average Eddington 
ratio for the active galaxies A!^ ™*^, which is still uncer- 
tain. IMcLure fc Dunlopl (|2004[ ) estimated that the av- 
erage accretion rate m|™'' varies from 0.1 at z ^ 0.2 
to 0.4 at z ^ 2 from a large sample of SDSS quasars. 
iKollmeier et al.l (|2006[ ) estimated the black hole masses 
and Eddington ratios for a sample of luminous AGNs 
with 0.3 < z < 4, and found that their average Edding- 
ton ratio is ~ 0.25. We adopt \^^ — 0.25 in this work in 
the calculation of active black hole mass density p^Yi{z) 
and then the inactive black hole mass density p™-^^'^{z)- 
The derived active black hole mass density p'^^{z) is pro- 
portional to 1/A|™'' (see Eq. [9]). We find that the derived 
pg^'(z) < piTH-z) for z < 3 (see Fig. [2]). The inactive 
black hole mass density /3{,'Sf'^*(z) is calculated by sub- 
tracting the active black hole mass density Py*!-^) from 
the total black hole mass density pbh(-2) at redshift z (see 
Eq. [TU|) . and /9bh(^) is derived with QLF as required to 
match the measured local black hole mass density pIJ^^^' 
at z = 0, which implies that the derived p™-^'^^{z) is al- 
most insensitive to the value of A^c^ for -^ ^ 3. The 
XRB and growth of massive black holes are mainly con- 
tributed by the accretion in the sources at low redshifts, 
which implies that our main conclusions on the growth of 
massive black holes at their late stage will not be affected 
if the value of A^ct ' does not deviate much from the value 
adopted in this work. For the average radiative efficiency 
e = 0.134 derived in this work, the average mass accre- 
tion rate for active galaxies rrV^^ = 0.187 corresponding 
to A|™'' = 0.25. The black hole mass densities can be 
calculated alternatively by using the Eddington ratio dis- 
tribution given by Eq. 1121 This distribution shows that 
most active galaxies have mass accretion rates m close to 
the critical one merit = 0.01 with an average m^^ ~ 0.1 
for Km = 0.3. As discussed above, our main conclusions 
will not be altered either with this value or a distribution 
instead of a single average one. The contribution of all 
bright AGNs to the XRB can be directly calculated with 
the QLF, w hich is independent o f the Eddington ratio 
distribution. IHopkins et al.l (|2007f) 's calculation showed 
that the observed XRB can be roughly reproduced with 
their QLF. 

Our calculations of the XRB contributed by ADAFs in 
inactive galaxies are sensitive to the spectra of ADAFs, 
which are described by several parameters, i.e., magnetic 
field strength relative to gas pressure /?, the fraction of 
energy directly heatin g the electr ons (5, and the viscos- 
ity parameter a. In I Caol ()20Q7D 's work, the theoreti- 
cal/observational constraints on the values of these pa- 
rameters were discussed, which leads to narrow ranges for 
the values of these parameters. The X-ray/7-ray spec- 
tra of ADAFs are almost independent of the value of /3 
(magnetic field strength), which only affects the spectra 
in radio bands. As we are focusing on the hard X-ray and 
7-ray energy bands, we find that the spectra of ADAFs 
depend most sensitively on the value of (5, because the 
electron temperature in the ADAF is sensitively affected 
by the energy directly heating the electrons. Besides a 
conventional value of ^ = 0.1 adopted in our calculations, 
we also carry out the calculations with 8 — 0.01 for com- 
parison. In this case, the radiative efficiency of ADAFs 



is lower than that for (5 = 0.1, as the heating of electrons 
is significantly suppressed (see Figs. |4]|6l). 

Figure [6] shows how the Eddington ratios -Lboi/^Edd 
vary with mass accretion rate m for different values of E 
and black hole spin parameter a. We find that the Ed- 
dington ratio varies with m roughly as iboi/-^Edd cc 'ti'*, 
where s ~ 1.2 — 2.1 depending on the values of parame- 
ters adopted. For comparison, we plot several power-law 
lines iboi/-^Edd oc rrf ^ s = 1, i.e., constant radiative effi- 
ciency e = 0.1 (black dash-dotted line), s = 1.2 and 2.1 in 
the same figure. The result s ~ 1.2 — 2.1 means that the 
radiative efficiencies of ADAFs increases with m. Our 
results are consistent with those adopted in the previous 
works (jNaravanfcYil 119951 lJesteHl2005t IHopkins et~all 
l2007f ). It is not surprising that the radiative efficiency 
increases with 5, and the bolometric luminosity Lboi is 
higher for rapidly spinning black holes provided the same 
accretion rate m is adopted. 

Combining the derived inactive black hole mass den- 
sity Pbif'^'(z) and the average spectra of a population 
of ADAFs with the mass accretion rate distribution 
(Eq. [12]), the contribution of the ADAFs in all inac- 
tive galaxies to the XRB can be calculated. The re- 
sults are plotted in Figs. [5] and [SI which are com- 
pared with the observed XRB. The contribution of bright 
AGN s to the XRB was estimated in many previous works 
(e.g lUcda et al.|[2003t iTreister fc UrF;^f200a iGiiTi et al.l 
[2003 LTrcistcr et al.l I2009H " We compare the sum of 
the contributi on of bright AGNs to XRB estimated by 
iTreister et"all ({2009 ) and that of the ADAFs in all in- 
active galaxies calculated in this work with the observed 
XRB. It is found that n^ > 0.3 for d = 0.1 (or k^ > 0.2 
for 6 = 0.01) is required in order to let the calculated 
XRB not surpass the observed XRB (see Figs. [5| and 
[9|). The dependence of required value of k^ on param- 
eter 6 can be understood that, a smaller 6 corresponds 
to a lower radiative efficiency, while a smaller Km repre- 
sents a flatter dimensionless mass accretion rate distribu- 
tion which has relatively more high- m(sa 10~^ — 10~^) 
sources. The average mass accretion rate for inactive 
black holes is relatively high for a distribution with a 
small Km, which also corresponds to a relatively high av- 
erage radiative efficiency, because the radiative efficiency 
of an ADAF increases with m (see Fig. ^. This implies 
that the distribution with a too flat mass accretion dis- 
tribution (i.e., small Km) will over-produce sources ac- 
creting at rates close to rhmax, which will radiates too 
much to the XRB. On the other hand, the value of Km 
can also be constrained by k and s according to the 
relation Lboi/^Edd oc mfand equations pTjl and P^ . 
IHopkins fc Hcr nguistl (|2009f ) suggested k « 0.3 — 0.8, and 
our calculation predicts s ~ 1.2 — 2.1 for different param- 
eters, then we can expect Km ~ ks ^ 0.36 — 1.67, which is 
roughly consistent with the constraints from comparison 
with the XRB in this work. Our results implies that the 
XRB contributed by inactive black holes is dominated 
by the radiation from high-rh. sources, rather than low-m 
sources. 

The energy peak of the contribution of the ADAFs in 
inactive galaxies to the XRB is around ^ 100 — 200 keV 
depending on the values of i5, which accounts for ~ 15 — 
20% of the XRB at these energy peaks (see Figs. [5| and 
[9]) . The Swift measurements on the XRB have estimated 
errors of ~ 3% (.AicUo et al..,200 8; .Trcister e t al..,2009.) . 
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Thus, we suggest that more accurate measurements of 
the XRB in the energy band with > 100 keV in the 
future may constrain the growth of massive black holes 
at their late stage more precisely. We note that the peak 
energy of the contribution of Compton-thick AGNs is 
~ 30 keV, while it is - 100 - 200 keV for that of ADAFs 
in inactive massive black holes. This implies that the 
constraints on the growth of massive black holes at their 
late stage from the XRB have hardly been affected by 
the possible uncertainty of the space number density of 
Compton-thick AGNs. 

The values of Kni=0.2 and 0.3 correspond to aver- 
age mass accretion rates of inactive galaxies rhf™^^. = 
8.35 X lO^'^ and 6.12 x lO"'', respectively (see Fig. H]). 
This means that the mass fraction of local black holes 
grown in ADAF phases should be less than 10.9% (for 
S = 0.1) and 14.5% (for S = 0.01) (see Fig. HJ. Our 
res ult i s abou t twice of that (^ 5%) for S = 0.1 given 
in ICaol (|2007l ). The discrepancy is mainly attributed 
to two factors. The first one is th a t a si ngle average 
mass accretion rate adopted in iCaol ()2007D for the cal- 
culations of ADAF spectra and their contribution to 
the XRB, which over-estimated the average radiative ef- 
ficiency for ADAFs in inactive galaxies, and therefore 
under-estimate the fraction of local black hole mass ac- 
creted during ADAF phases. The another one is that a 
different AGN population synthesis model for the newly 
measured XRB with Swift is adopted in the calculations 
in this work. Our present results are roughly consistent 
with that derived from observed Eddington ratio distri- 



butions in lHopkins et al.l (|2006D . 

ADAFs are very hot, and the temperature of the 
ions in ADAFs can be as high as ~ lO^^K (see Fig. 
[3|), which implies that 7-ray emission may be pro- 
duced through the pion production processes in the 
proton-proton (p- p) collisions and subseq uently decay 
of neutral pions (JMahadevan et al.l I1997D . We calcu- 
late their contribution to the extragalactic 7-ray back- 
ground (EGRB), and find that less than 1% of the ob- 
served EGRB is contributed by the ADAFs in these faint 
sources if the massive black holes are spinning rapidly, 
while the contribution to the EGRB from ADAFs can 
be neglected if the black holes are non-rotating (see 
Figs. [To] and [TT|) . Our results are consistent with 
the previous works showing that about ~ 25% to ~ 
100% of the EGR B can be attributed to the unre- 
solvcd blazars (e.g., Padovani ct al. 1993; Chiang ct al 
1995: Stecker fc Salamon 1996; Miicke fc Pohl :200' 
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